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/\ ' Abstract 



We study the end-to-end distribution function for dilute polymers. We 
present a computation to order O(e^), e = 4 — d, and discuss in detail its 
asymptotic behaviour for small and large distances. The theoretical predictions 
are compared with Monte Carlo results, finding good agreement. We show 
that the McKenzie-Moore-des Cloizeaux phenomelogical ansatz provides a very 
precise approximation to the exact EEDF. 



1 Introduction 

The statistical properties of dilute polymers in good solvents have been the subject of 
extensive studies during the years [|1]-§1 • A significant understanding of the problem 
was reached when it was realized that long polymers could be modelled by chains 
with an excluded-volume interaction. This allowed the introduction of simplified 
theoretical models which could be analyzed more easily. From a theoretical point 
of view, an important step forward was made by de Gennes |^, who proved that 
the statistical properties of polymers could be obtained as the limit A^ — > of the 
iV-component 0^ theory, opening the field to the many methods that were developed 
at the time for the study of the critical behaviour of spin systems. 

In this paper we consider the end-to-end distribution function (EEDF). During the 
last thirty years a lot of work has been devoted to the study of this quantity. Exact 
results were obtained in Refs. [|6|-[TT[], and many numerical studies checked some of 



the theoretical predictions p!^|-[2D|]. In particular, these works tried to understand 
which of the several phenomenological expressions pl| provided the best description 



of the numerical data. The proposal by Mazur ||^ was clearly excluded ||T5[, while 
the theoretically motivated form of McKenzie and Moore and des Cloizeaux 0, 
TI| was confirmed to a good accuracy |T^, |l^, |20| . It is interesting to notice that a 



precise knowledge of the EEDF could be relevant in experimental studies. Indeed, 
as observed by des Cloizeaux and Jannink 0, the EEDF could be determined from 
scattering experiments with a dilute solution of polymers with marked endpoints. A 
measurement of the scattered intensity at large angles would determine the EEDF in 
the large-momentum region. This would provide an estimate of the critical exponent 
7, which is otherwise inaccessible experimentally. 

In this paper we reconsider the problem of the determination of the EEDF. We 
extend the e-expansion calculations of Ref. |23] to order e^ and give perturbative 



expansions for the quantities that characterize the asymptotic behaviour for small 
and large distances. Using the e-expansion we can show that the phenomenological 
parametrization of |P,|lll is essentially exact for distances much larger than the cor- 
relation length — the expected discrepancy is of order 1-5% in the region accessible 
to simulations — while in the opposite range the discrepancy should be (at most) of 
order 10%. We also give theoretical expressions for several quantities derived from 
the EEDF. Using precise Monte Carlo estimates of the critical exponents, we derive 
accurate theoretical predictions for the asymptotic behaviour of the EEDF. The the- 
oretical estimates are compared with numerical results obtained from a simulation of 
self- avoiding walks on a cubic lattice, finding good agreement. 

The paper is organized as follows: in Sec. 2 we introduce our notations and 
definitions and review the exact results that are available for the EEDF. In Sec. 3 
we report our computation of the EEDF to order e^. Only the results are given, the 
technical details being presented in the Appendix. In Sec. 4 we obtain estimates 
for the asymptotic behaviour of the EEDF using the Laplace-de Gennes transform. 
Finally in Sec. 5 we discuss the numerical results. 



2 Definitions 



We consider a inonodisperse ensemble of polymers with A^ monomers. If r is the 
vector joining the endpoints of the walk, we will be interested in computing the 
unnormalized distribution CAr(r) [^] of the endpoint vector. We also introduce a 
normalized distribution 

P«(r) = #i<4- , (2.1) 



ErC7v(r) 

the mean squared end-to-end distance 



and the related correlation length 
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^N — nJ^e,N- 



(2.2) 



(2.3) 



In the limit A^ — i> oo, |r| ^ oo, with |r|A^ '^ ~ |r|/^7v fixed, the function -PAr(r) has 
the scaling form [^,||,|lT 






(2.4) 



where p = r/C,N, P = \p\, d is the space dimensionality, and A is a correction-to- 
scaling exponent. By definition 



p'-'dpfip) 
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where Sd is the volume of the d-dimensional sphere 

27r'^/2 

Od 



(2.5) 
(2.6) 

(2.7) 



Tid/2) ■ 
Several facts are known about /(p). For large values of p it behaves as P,r7|,P|,pT| 

/(p) ^/ooP'^exp(-Dp^) , (2.8) 

where a and 6 are given by 
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For p — * 0, we have P,P 



1-z/' 
2ud-2-f + 2-d 

2(1-^) 



f(p) ^ /o 



(2.9) 
(2.10) 
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where 

e= Tzl _ (2.12) 

We can also consider the Fourier transform of / (p) , 

which is the critical (large- iV) limit of CAr(q)/cjv(0) with Q = q,^jv fixed, CAf(q) being 
the Fourier transform of cjv(r). For Q — *^ 0, f{Q) has a regular expansion in powers 
of Q^, while for Q^ — ^ oo it behaves as 

KQ) = fooQ-'-'. (2.14) 

The constants foo and /o are related. Indeed 

/o^Z-C")-"'^. (2-15) 

A phenomenological representation for the function /(p) has been proposed by McKen- 
zie and Moore and des Cloizeaux 0: 

/(p) ^ /ph(p) = Up' exp (-Dphp') . (2.16) 



Here 6 and ^ are fixed by ( p.9| ) and ( |2.12| ), while /ph and Dph are fixed by the nor- 
malization conditions (G.5|) and (|2l 



In two dimensions 



_ r[(i-z.)(2 + rf + g)] p 
^'^^ - i2dr[(i-.)(d + ^)]/ ' ^2-^^^ 



(2.19) 
(2.20) 





5dr[(i-z/)(d + e)] 


l^ph 


= 0.026339478... 


/ph 


= 0.046757638... 



while in three dimensions, using the precise estimate u = 0.58758 ± 0.00007 ||2^ and 
our result [^] 7 = 1.1575 ± 0.0006, we have 

Dph = 0.14470 ± 0.00014, (2.21) 

/ph = 0.015990 ± 0.000008, (2.22) 

9 = 0.2680 ±0.0011, (2.23) 

6 = 2.4247 ±0.0004. (2.24) 

Notice that /ph(p) cannot be exact since the exponents 6 and a are different, but 
this is less crucial in rf = 3, since their numerical values are very similar. In order 



to improve the approximation, one can also take 9 and 5 as free parameters, and 
determine them by fitting numerical estimates of /(p). However, once 9 and 5 are 
given, Z)ph and /ph are uniquely determined by ( p.l7| ) and ( p.l8| ). In the following we 



will indicate with "phenomenological representation" the function (|2.16|) with 5 and 
9 fixed to their theoretical values. 

For the purpose of computing D and 5 from Monte Carlo data, it is much easier 
to consider the "wall-to-wall" distribution function 

Pn,w{x)= Y1 PN{x,X2,...,Xd), (2.25) 

X2,...,xa 

which represents the probability that the endpoint of the walk lies on a plane at a 
distance x from the origin of the walk. In the large- A^ limit, Pj\ry^{x) has the scaling 
form 

PnA^) = -^fM (l + 0(iV-^)) , (2.26) 

where p = x/C,n- We will show in Sec. § that for large p we have 

fM - UooP"- exp(-Dp^) , (2.27) 

where 5 is given by (p.9|), D is the same constant appearing in Eq. ( p.8|) , and 



(y^ = 5{v-^ + ]^. (2.28) 

3 Distribution function to order O(e^) 

We will now derive the EEDF /(p) using a continuum description and the standard 
e-expansion. We start from the Edwards' path integral [^ for the probability distri- 
bution function of the end-to-end distance r of a chain with contour length N m. d 
space dimensions 

^W=^ _ f 1 /-^ , (dMs)V w r^ , r^ 



(3.1) 

Here x(s) is the position vector of the arc-length position s of the chain, the integral 
over D\x\ represents the summation over all possible configurations between the two 
ends of the chain, and w is the unrenormalized strength of the excluded-volume 
interaction. We use adimensional units setting the Kuhn step length equal to 2d. 
In this Section we report the computation of the EEDF /(p) to order e^, where 



as usual e = 4 — d, extending the results of Ref . p2| . The diagrams that need to be 
computed are reported in Fig. |^. We obtain: 

• Graph (a) 



1 „ /„ e _ t 



2 V 2 2 
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(3.2) 



Graph (b) 



• Graph (c 



4 V 2/ V 2' 2 

-Nq 



X 



l-Nq')e~'"^ +2eg,{Nq')+Ae'g2{Nq') +0(e); (3.3) 



1 / e e 

-5 2- -,-l + - 
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iv d ' 



1 _ Nq')e-'''^' + 2e^i(iVg^) + 4e^^2(iVg') + 0(e); (3.4) 



• Graph (d) 
Gd{Nq^) = 



3 3 3 7r2 r2 3, 3, o ' 

— 1 loer 2r H — loer 2r 

2e2 4e 2 8 6 2 ^ ^ 2 ^ 



X 



1 - Arg2)e-^'?' + 2e^i(Arg2) + 4e2^2(A^g') 
+AN2{Nq^) + 4Ar3(Arg2) + o(e). (3.5) 



The functions gi{Nq^), g2{Nq^), g^iNq^), N2{Nq'^), and N^iNq"^), are defined in the 
Appendix. The result for diagram (d) depends apparently on an arbitrary parameter 
r; note that the functions N2{Nq'^) and N^lNq"^) depend also on r in such a way to 
make the final result independent of r. In principle r can be set to any value. We 
have kept it arbitrary, in order to have a check of the calculations: indeed the final 
results must not depend on r. At two loops we obtain for the unrenormalized CAr(q) 



Cat! 



q) = e-^^' + wN'/^GaiNq') + w^N' \Gb{Nq^) + G,{Nq^) + Gd{Nq^)\ + 0{ 



[w 



(3.6) 

where w = wN^ and A^^^ = 2(47r) '^/2/r((i/2). The computation of the universal 
EEDF goes through several steps. First of all, we compute the correlation length ^^v 



using the definition (|2.3| ) and then we express A^ in terms of S,n- We obtain 



A^ = ^^ l + aiw Cn + «2W^ i^j^ + 0{w 



where 
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+ — log 2r + :^ log 2r - - log^ 2r + -(A^2(0) + A^2(0))- 
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(3.7) 



(3.8) 



(3.9) 



where 7^ ~ 0.5772156649 is the Euler's constant. 

Substituting this expression into Eq. (|3.6|) , C7v(q) becomes a series in w^%, 
with coefficients depending on the combination q^^v = Q- One can then compute 
f{Q) = CAr(q)/c7v(0). This quantity, once expressed in terms of Q, requires only a 
renormahzation of the interaction strength w in order to be finite. In the minimal 



subtraction scheme we have |3,27| 



W = Wr[1 + -Wr + 0{Wr^ 



(3.10) 



The expansion in terms of wr is finite. The critical theory is obtained replacing wr 
with its fixed-point value w*r, 



The final result can be written in the form 

KQ) = e~Q' + ep\Q) + e'p\Q) + O(e^) 
where 



(3.11) 



(3.12) 



(3.13) 



P\Q) = {^-l^)9iiQ') + ^UQ')-,l;sUQ') + 7mQ') + WQ') 
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128 32, 

1 7i?\^2w, 



16' 



128' 



TC 



2\2-Q^ 



768 



m'e 



,16 si)^'^'^^^'^^^^'^^^' 
where, for any function h{Q'^), we define 

h{Q^) = hiQ') - e-«'/i(0)(l + Q') - e-'^"h'{0). 



(3.14) 



(3.15) 



From the results reported in the Appendix one can derive the asymptotic behaviour 
of f{Q) in the limits Q — > and Q -^ 00. For Q ^ we have 



m 



i-(Q 



2\2 



where we have introduced 62 

&2 



2 + 32+ 1024 + 768 + Yh+^(^ 



+ 0{Q' 



dxx^Kf{x) ^ -0.000235007, 



1536 Jo 
where Ki{x) is a Bessel function. For Q^ — > 00 we have 



(3.16) 
(3.17) 



m 



1 



{Q 



2^2 



.l + |,(ll-12,,)-|!logQ^ 



+ 0{e' 



(3.18) 



where terms of order log Q"^ / {Q'^y have been discarded. The function /(p) can be 
derived from f{Q). We will be interested in its asymptotic behaviour for p -^ oo. 
Using the results of the Appendix, a lengthy computation gives 



/(P) 



1 e 

+ 



167r2 647r2 



5127r2 



log^ + 21og47r+|7£;-l) 



log' ^ + (41og47r + 7^ + i) log ^ + 1 + 



i9vr2 
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+s7s + 171-4 log 47r + 2^e log 47r + 4 log' 47r + A A 



(3.19) 



where terms vanishing for p —>■ have been neglected. The constant A is defined 
in Eq. ( A.43|) . Numerically A ^ —1.30204. The expressions ( |3.18 ) and ( 3.19|) are 
compatible with the asymptotic behaviours ( p^.l4D and ( |2.11D and give the following 
expansions for the constants foo and /q: 



/oo = -|e+2l6e'(ll-127ir) + 0(e3) 



(3.20) 



/o = T?^ + ^(21°g4vr + i7.-l)+34;,(l + ivr' + i7. + hl 






-4 log 47r + 27£; log 47r + 4 log' in + iA) + 0{e 



(3.21) 

A check of these results is provided by Eq. ( |2.15| ). It is easy to verify that the two 
previous expansions satisfy Eq. ( |2.15D to order e. Moreover, using Eq. ( |2.15| ), we can 
obtain [EH] the contribution of order e^ to foo- Explicitly we have 



(- 



41 I 21 
'4096 ~'~ 1024 '^ 



1024 



1e 



2048 



TC 



+ 1C(3) - hA) e\ 



(3.22) 



Since only three terms are available, it is difficult to obtain reliable estimates from 
these expansions. Setting simply e = 1, we obtain 



/o 



-0.0557, 



/o ^ 0.0176. 



(3.23) 



which takes into account the 



We have also used the resummation method of Ref. 

singularity structure of the Borel transform of the perturbative series. We obtain |^| 

/^ ?« -0.089 ± 0.017, /o^ 0.015 ±0.002. (3.24) 

The ratio /oo/Zo can be determined precisely using Eq. ( |2.15| ) and the precise Monte 
Carlo determinations of the critical exponents: 



f 

^ = -4.883 ±0.017. 

Jo 



(3.25) 



We can use this result to test the accuracy of the estimates ( p.23| ) and (|3.24|) . Using 
Eq. ( p:23D we obtain f^o/ fo = -3.16, while Eq. (^2§ gives foc/fo = -5.9 ± 1.4. 
The estimates have an error of approximately 20-30%. It is reassuring that the errors 
obtained with the resummation method correctly describe the discrepancy. 



4 Asymptotic behaviour via Laplace-de Gennes trans- 
form 

Improved estimates of the asymptotic behaviour of /(p) can be obtained from the 
precise results that have been derived for 0{N) spin models, using the Laplace-de 
Gennes transform method p,|^, |lT|. We introduce in d dimensions the two-point 



function 

oo 

G{r;(3)= Ee-^^c^(r), (4.1) 

N=0 

its Fourier transform 

G{c^;(3)=J2e''^-^Gir;(3), (4.2) 



and the correlation length 



^ ^^^- 2d E.G(r;/?) ^^-"^ 



In the critical limit G(r; (3) and (^(q; (3) have the scaling form 

«(-« .. 1 ^Oi,,), M«5(Q). (4.4, 



T.rG(v;P) £,'(13) '"' G(0;/3) 

where p = r /^{(3) and Q = \q\^{P). The function D{Q) is the d-dimensional Fourier 
transform of D (p) . 

Several facts are known about the function D{Q). For Q ^ 1, D{Q) has a regular 
expansion in terms of Q^ as 

oo 
j=2 



The coefficients 6„, n = 2,3,4,5 have been computed [Q in the e-expansion up to 
O(e^) and [Q in the fixed-dimension expansion in d = 3 up to 0{g^). It turns out 
they are extremely small and satisfy &2 ^ &3 ^ &4- • •• For 62 the explicit formulae 
are 

62 = -0.000235007e2[l + 1.0632e + O(e2)], (4.6) 

62 = -0.00015432/[1 + 0.0780213^ + 0.0465896/ + 0(/)]. (4.7) 



Here g is the renormalized four-point coupling |^ constant whose critical value pO 
37|-^, in the normalization we use, is g* ~ 1.39. The expansions for 62, &3, • 



can be resummed using the method of Ref. |^, obtaining |^ 62 ~ (—3 ± 1) ■ lO"''. 
Analogously 63 ~ 10~^. An exact-enumeration study confirmed these results and 
gave the bound —3 ■ 10""^ ^ &2 ~ 0. In two dimensions estimates have been obtained 
from the analysis |^ of exact-enumeration expansions on the triangular, square and 



honeycomb lattices: 62 ~ 0.00015(20), I63I ^ 3 ■ 10 ^. In the following we will not 
need the explicit values of the constants 6j. Indeed they are too small to give any 

9 



numerically important effect. Thus, for Q ^ 0, we can approximate D ^{Q) with 
1 + Q2. 

For Q ^ 1, the behaviour oi D{Q) is predicted by a short-distance renormalization- 
group analysis ||4^,0 and one has 



D{Q) = ^ + 



Di 



+ 



D, 



+ . 



Q2-»? Q2-n+{^-(i)/u Q2-n^\lv 

The exponents r] and a are related to the exponents 7 and v by 



7/ = 


= 2-2. 

V 


o; = 


= 2-du. 



1 - 0.03173916^ - 0.0353978e^ + O(e^) 
i (1 + ^6 -0.482262)+Di36^ + 0(6^), 

-^ (1 + ^e - 0.176536^) + D^^e' + O(e^). 



The constants Do; -Di and D2 have been computed in the e-expansion [p^ ,|34 

Do 
Di 

D2 

The coefficients Di^ and D23 are not known but satisfy the relation: 

L)i3 + D2Z = 0.05934. 
Resumming the perturbative series, we obtain 

Do -0.97 ±0.02, Di ^ 0.71 ± 0.04, D2 ~ -1-70 ± 0.06 



(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 

(4.13) 

(4.14) 
(4.15) 



Using the longer series for Di + D2 we obtain Di + D2 ^ —0.97 ± 0.02. It has been 



remarked by Aharony and Fisher |^3[ that one can also rewrite 

A = - (7 - 1 + 0{e 



(4.16) 



and thus, for d = 3, using the known values of the critical exponents, we obtain a 
similar estimate D^ ^ 0.67. 

Finally one can determine the large-p behaviour of D{p). For p ^ 00, using the 
notations of Ref . [Q , we have 



D{p) ^ A*p-^'^-''^/\- 



•S'lfP 



(4.17) 



where A* and S^j are non-perturbative constants. The constant A* can be related to 
the residue S^ of the propagator at the mass pole. It is given by 



St 



2 ag^ 



D-'{Q) 



(4.18) 



?2=-5* 
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One obtains 

A* = ^(27r)(i-^)/2 iSljf-'^^' . (4.19) 

The constants S*^ and S\,i have been computed [^] in the e-expansion up to O(e^) 
and in the expansion in fixed dimension d = 3 up to 0{^g^\. 

S*2 = 1 + 0.0004891006^ + 0.0005224256^ + 0(6^), (4.20) 

SIj = 1-0.0002412746^-0.0002573036^ + 0(6^), (4.21) 

S*z = 1 + 0.000335466/ + 0.0000140322/ + 0(c/^), (4.22) 

SIj = 1- 0.000163057/- 0.0000088348/ + 0(/). (4.23) 

It is evident from these expansions that both constants are one with very small 
corrections. Resumming the expansions, we obtain in three dimensions Sl^ ~ 1 — (3 + 
1) ■ 10~^ and S"^ ~ 1 + (5 ± 1) ■ 10"''. In two dimensions estimates have been obtained 
from the analysis of exact-enumeration expansions ^1[|: Slj ~ 1 + (1 ± 2) • 10~^, 
5^ ^ 1 - (2 ± 4) ■ 10-^ 

From the asymptotic behaviour of D{p) we obtain corresponding predictions for 
/(p), using the fact that P,p!T[ 



1 f+ioo 

CN{r)=— d(3e''''G{T;(3). (4.24) 

zm J-ioo 

We report here the results; the derivations can be found e.g. in P,pT|. 

We begin by computing the large-Q behaviour of f{Q). Using Eq. ( |4.8|) we obtain 
the expression (|2.14|) with 



~ 1 / r(7) \ ('+")/' 

f^ = -r(7)Z^i sm{7idu)T{du) L^^^i^u) ) ' ^^'^^^ 

Using the numerical values of the exponents, we obtain in two dimensions 

e = —, (4.26) 

f 

^ = -0.11062768... (4.27) 



and in three dimensions 



9 = 0.2680 + 0.0011, (4.28) 

^ = -0.12393 + 0.00026. (4.29) 



Using Eq. ( |2.15| ), we obtain correspondingly 



^^r0.050548 forrf = 2, 

Di~{ 0.02539 + 0.00014 for d = 3. 
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(4.30) 



In two dimensions we do not have any prediction for Di, but, on the basis of the e- 
expansion result, we expect Di to be of order 1, so that a reasonable guess is —0.05 < 
/oo ~ —0.15 and 0.02 ^ /o ^ 0.07. In three dimensions, using Di ^ 0.71 ± 0.04, we 
obtain /oo ~ —0.088 ± 0.006 and /o ~ 0.018 ± 0.001. These estimates are in good 
agreement with the results obtained in Sec. |^. 

It is interesting to compare the small-p behaviour of /ph(p) with that of the exact 
function /(p). Using Eq. ( |2.16| ), we obtain /o^ph ~ /ph2^ ~ 0.0193, which does not 
differ significantly from the estimate of /o reported above. We can also compare the 
two predictions within the e-expansion. We have 

-^ = 1 + O.109663e2 + O(e^), (4.31) 

which shows that the phenomenological approximation is essentially correct, with an 
expected discrepancy of order 10%. 

Using Eq. (|4.17|) , one obtains the asymptotic behaviour (|2.8|) where /oo and D 
are given by 

D = l^(.''S'„,,m^"'' (4.32) 



r(7 + 2z/ 



2S*z (27r)^/2(l - i/)i/2 V ^m) l^r(7 + 2z/) 



Using the values of Sl,i and S*^ we have reported before and the values of the exponents 
7 and v we get in two dimensions: 



5 = 4, 


(4.34) 


5 

^ = 8' 


(4.35) 


D = 0.02771 ±0.00001, 


(4.36) 


/oo = 0.04273 ± 0.00002. 


(4.37) 


In three dimensions, using u = 0.58758 ±0.00007 p4| and our result |25| 7 = 


= 1.1575± 


0.0006, we obtain 




5 = 2.4247 ± 0.0004, 


(4.38) 


a = 0.255 ±0.002, 


(4.39) 


D = 0.1434 ±0.0002, 


(4.40) 


/oo = 0.01581 ±0.00002. 


(4.41) 



Notice that the estimates of D and /oo would not have significantly changed, had we 
used the Gaussian values 5*2/ = 5*^ = 1. In three dimensions the error is dominated 
by the error on 7 and u. 

We can also use Eqs. ( [4.32|) and ( |4.33| ) to derive e- and (^-expansions for D and 



/oo- We obtain in 4 — e dimensions 

D = ^-0.0877837e-0.0327168e2 + 0.0394476e3±O(e^), (4.42) 

/oo = --r^±0.00579036e±0.00223513e2 + 0.00152100e=^±O(e^); (4.43) 
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in fixed dimension d = 3 we get 
1 



D = -- 0.0877837^ + 0.0153578^^- 0.010465V + 0(/), (4.44) 



/oo = —4^ -0.00534423c/ + 0.000981669/ -0.000734351c/^ + O(c/^). (4.45) 
(27r)'^/^ 



Resumming the expansions using the method of Ref. ^^, we obtain in the two cases: 
the e-expansion gives D = 0.1461±0.0040, /oo = 0.0138±0.006, while the (/-expansion 
gives D = 0.1445 ± 0.0015, /oo = 0.01602 ± 0.00010. These results are less accurate 
than the previous ones, that were obtained using the precise Monte Carlo estimates of 
the critical exponents. Note that the estimate of /oo obtained using the e-expansion 



is not compatible, with the quoted error bars, with the estimate ( [4.41|) . This is not 



surprising since in the expansion ( [4. 431) all coefficients have the same sign: therefore. 



a Borel resummation based on the large-order behaviour of the coefficients (that 
predicts coefficients alternating in sign) is not expected to work well. 

It is interesting to compare the estimates (|4.40| ) and ( 4.41 ) with the phenomeno- 



logical approximation ( |2.21D and ( |2.22| ). It is remarkable that the discrepancy is so 



tiny, precisely of 0.8% for D and of 1.1% for /oo. Also in two dimensions the phe- 
nomenological approximation works reasonably well: the discrepancy is of 5% for D 
and of 9% for /oo. This nice agreement can be understood within the e-expansion. 
Indeed 



D 

Too 
/ph 



1 - 0.012109e2 + 0.0039898e^ + O(e^), (4.46) 

1 - 0.0083917e2 - 0.0065622e3 + O(e^). (4.47) 



The e-expansions of D and /oo and of their phenomenological approximations Dph 
and /ph differ by terms that are very small. Setting e = 1, one finds an expected 
difference of order 1%, in agreement with the estimates above. 

Notice that also the exponent a does not differ significantly from 6 in three di- 
mensions. This explains the success of the phenomenological approximation ( p.l6| ) 
for p -^ oo. Indeed, in this limit, we have 



/ph(p) ^ /ph^.-.^^p 



f{p) /c 

so that 



{Dpi,-D)p' ^ 1.011 p""''exp(-0.0013p'), (4.48 



fjp) - fAp) 

fip) 




for p = 2, 

for p = 5, (4.49) 

for p = 10. 

However /(5) ~ 2 -10"^, and /(lO) ~ 8-10^^^, so that in practice /(p) can be sampled 
up to p ~ 5-6. Therefore, in the region accessible to numerical simulations, /ph(p) 
provides an accurate description of the large-p behaviour of the EEDF. 
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Finally, we can use ( [4 .51) to get predictions for the moments of /(p). It is simple 
to show that the invariant ratios 



Er \r\'"'PN{r) 
{Er \r\'PNir)) 
approach, for N ^ oo, universal constants M2f^ given by 



M.. = .t:':i:"::L (4.80) 



^^- - r(.''l.t)rw^- 1^ - ''■^<* - ^)1 *' n (1 + 1) . (4.51) 

where, we have neglected all bn with n > 3 and all powers of 62- Notice that Eq. 
( [4.51| ) is exact for k = 2. 



Again, we can compare the exact expression ( f4.51| ) with the prediction obtained 
by using the phenomenological representation ( |2.16| ): 



Mi... = r ^^^j r (^— j r^^— j . (4.52) 

Numerically we have in three dimensions: 

Ml = 1.51397(79) M*^^^ = 1.50876(23), (4.53) 

M* = 3.018(4) Mp^e = 2.993(1), (4.54) 

M* = 7.392(15) Mp^s = 7.292(5), (4.55) 

M*Q = 21.35(6) Mp\ 10 = 20.94(2). (4.56) 

The phenomenological predictions are definitely not exact, but they show very small 
differences with respect to the exact ones. Since the moments define uniquely the 
distribution function, this implies that Eq. ( |2.16|) is not only a good approximation 



in the regions of large and small values of p, but that it also provides a good overall 
parametrization of the EEDF. 

We can also consider the "wall-to-wall" distribution P^^ix). It is easy to see that 
for large p we have Eq. ( p^.27D with 



xf 
a^i = u - 

V 


-^+D' 


Jw,oo loo {■^ 


' [ni+2u) 


In two dimensions we obtain 




CTw 


3 

8' 


Jw,oo 


= 0.32167 ±0.00017 


and in three dimensions 




Cw 


= -0.1695 ±0.0016, 


Jul, 00 


= 0.2855 ± 0.0004. 



5{l~d)/4 



(4.57) 
(4.58) 



(4.59) 
(4.60) 



(4.61) 
(4.62) 
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5 Monte Carlo study of the distribution function 

The EEDF has been extensively studied numerically in three dimensions. The Monte 
Carlo work essentially focused on the exponents and verified that the data could be 



com- 



(5.1) 



well described by the phenomenological expression (|2.16|) . Baumgartner |T4 
puted the exponent S obtaining 

(5 = 2.44 ± 0.05, 

in very good agreement with the theoretical result ( |4.3^ ). The exponent 6 appearing 
in the phenomenological expression ( |2.16|) was computed by several groups, obtaining 

0.270 ±0.006 Ref. [14], 



0.27 Ref. 

0.262 ± 0.013 Ref. 
0.224 ± 0.006 Ref. 



0, 



m 



(5.2) 



m- 



All estimates but the last one do not differ sensibly from our theoretical results for a 
and 9, and as expected, they lie between these two estimates. The estimate of Ref. |jl9 



is instead too low; probably, the numerical data are affected by large corrections to 
scaling. 

In this Section we want to extend these numerical analyses, checking the renor- 
malizat ion-group predictions presented in the previous Sections. We will not use the 
phenomenological expression (|2.16|) , but we will compute the exponents 6, a, and 9 
studying the large-p and small-p behaviour of the EEDF. At the same time we will be 
able to compute the constants /o, /oo and D and to compare them with the theoretical 
predictions. 

In order to compute the EEDF, we have generated A^-step self- avoiding walks on 
a three dimensional cubic lattice, using the pivot algorithm Hl-E^. Since in three 



dimensions corrections to scaling are particularly strong, we generated long walks 
with 500 < A^ < 32000. 

First of all, we have checked the prediction ( 4.51|) for the invariant ratios M2k, 
for fc = 2, 3, 4, 5. In Table |I| we report the Monte Carlo estimates of M2k for various 
values of A^. We have performed the extrapolation to A^ — > oo using 

A 



Mr 



2k 



M* + B; 



2k 



^000 

"aT 



(5.3) 



The final estimates are compatible with the less precise results of Ref. [^ and are 
in very good agreement with the theoretical predictions. The corrections to scaling 
appear to be quite important: we estimate A ^ 0.55 ± 0.08, in agreement with the 
renormalization-group and Monte Carlo predictions [Q. Note that, at our level of 
precision, the phenomenological predictions ( [4.52|) are not consistent with the numer- 
ical data. 

It is interesting to observe that we can use Eq. 



m 



to obtain independent 
= -(3 ± 1) ■ 10-^ 



estimates of the critical exponents. For instance, using 62 = 

u = 0.58758 ± 0.00007 g], and the Monte Carlo result for M'*, we obtain 



7 = 1.1576 ±0.0013, 



(5.4) 
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which is in perfect agreement with the more precise estimate of Ref. |25|. It should 
be noticed that most of the error on 7 is due to the error on M4, and, indeed, with 
the existing estimates of v and 62, it would be possible to obtain with this method 
an estimate of 7 as precise as that given in Ref. [^. By considering M^ and Mg one 
can try to estimate v and 7 simultaneously. One obtains v = 0.588(6), 7 = 1.159(24). 
Finally, using the estimates of 7 and i> and our Monte Carlo result for M4 , we obtain 
a bound on 62- We get I62I < 1-4 ■ 10~^. 

In addition to M2k, one can consider the invariant ratios 



^- = t:':.^\^:.i:\i (5.5) 



The rotational invariance of the critical limit gives 



In Table ^ we report the estimates of the ratios R2k for k = 2,3,4,5 and several values 
of A^. Notice that, in this case, corrections to scaling are practically absent. This 
is in agreement with the analysis of Ref. ||5^ which showed that quantities like i?2fc 
have corrections of the form A^""^ with r ^ 2v. Thus, they are much smaller than the 
standard corrections which behave as A^~^ ~ A^""'^. 

Let us now consider the EEDF itself. The computation of P/v(x) from the Monte 
Carlo data is straightforward, but it is less clear how to estimate the error bars. If 
the walks are generated independently, and p(x) is the estimate of P/v(x), the error 
is 



-f>{x){l — p{x)) 



^ ..-..- ..-.. , (5.7) 

where Nmc is the number of Monte Carlo iterations. In our case, however, the 
walks are not generated independently. Therefore, one should take into account the 
autocorrelation time of the algorithm and the fact that estimates at different values of 
X are correlated. In practice it is not feasible to take into account all these effects. We 
have simply observed that since global observables decorrelate after a few accepted 
pivot moves, a reasonable estimate of the errors can be obtained by replacing in 
Eq. ( pr?!) Nmc with /nNmc, where /a? is the acceptance fraction of the algorithm. 
Correlations between different points are neglected. 

The functions /(p) and fwip) are reported in Figs. || and 0. The data fall on a 
single curve as expected: within the accuracy of the plot, no corrections to scaling 
are visible, but, as we shall discuss later, corrections are present if one looks at the 
data in more detail. 

Let us now study the asymptotic behaviour of the EEDF. We will begin by con- 
sidering the wall-to-wall EEDF fw{p). In order to study its large-p behaviour, we 
have performed two different sets of fits: 

log(/.(p)) = log U,^- Dp', (5.8) 

log(p°-''V^(p)) = log^,oo-V- (5.9) 
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In the first fit we have neglected the power term p"^ that appears in the asymptotic 
behaviour of fw{p), in the second one we have taken this term into account using the 
theoretical prediction for a^, Eq. ( |4.61| ). For p ^ oo, both fits should give the correct 
result for D and 6, while only the second one gives an estimate of fw,oo- There are two 
types of systematic errors in these fits. First, there are scaling corrections: the scaling 
curve is obtained only in the limit N ^ oo. Secondly, there are non-asymptotic 
corrections: Eqs. ( |5.8D and ( |5.9| ) are valid only asymptotically for p — »• cxo. In order 
to detect scaling corrections, we have performed fits at fixed values of A^. Then, we 
have compared the results, looking for systematic variations of the estimates with 
the length of the walk. The final result is obtained including in the fit only walks 
with A^ > iVmin, where A^min is chosen so that the estimates for all A^ > A'min are 
independent of A^ within error bars. A similar strategy has been used to detect non- 
asymptotic effects: we have performed several fits using in each case only data with 
P > Pmin- Looking at the variation of the estimates with Pmin, "we can estimate the 
non-asymptotic corrections. The results of the fits for fixed values of A^ and for three 
different values of pmin are reported in Table 0. Apparently they do not show any 
systematic dependence on A^, except perhaps A^ = 500 and A^ = 1000: indeed the 
estimate of 6 for A^ = 500, 1000 are slightly higher than the estimates obtained for 
larger values of A^, while the estimates of D are slightly smaller. One may suspect that 
these results are affected by scaling corrections of size comparable with the statistical 
error. For this reason, our final estimates are obtained using all data with A^ > 2000 
only. We obtain from the first fit 

^ r 2.413 ±0.006 Pmin = 3.0, 

12.420 ±0.014 Pmin = 3.5, ^ ^^ 

^ r 0.150 ±0.002 Pmin = 3.0, 

10.148 ±0.004 Pmin = 3.5. ^ ' 



The second fit gives 



2.458 ± 0.006 Pmin = 3.0, 

2.455 ±0.013 Pmin = 3.5, 



(5.12) 



^ r 0.136 ±0.002 Pmin = 3.0, 

1 0.137 ±0.004 Pmin = 3.5. ^^^ 

These results are in good agreement with the theoretical estimates reported above, cf. 
Eqs. ( |4.38|) and ( [4.40| ). Indeed, increasing Pmin, we observe the expected convergence 



to the theoretical results. The second fit gives also an estimate of /^,oo: we find 

^r 0.274 ±0.002 Pmin = 3.0, 

Jw,oc \ 0.276 ±0.006 Pmin = 3.5, ^ ' 



which converges to the estimate ( [4.62| ) for large values of Pmin- 

In order to estimate a^ and /^,oo, we have performed fits of the form 

log /^(p) exp{Dp^) = log fnj,oc + CTyj logp, (5.15) 

using the theoretical estimates of D and 6, for various values of p > pmin- These fits 
are extremely unstable. Indeed the EEDF drops rapidly to zero (see Fig. 2) so that 
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the fit uses data in a small interval in which log p does not vary significantly. Results 
with reasonable errors can be obtained only for pmjn < 3, and thus we have analyzed 
the data with pmin = 2, 2.5 and 3. The results for fixed values of N are reported 
in Table |[ Looking at the table, one immediately sees that there are strong non- 
asymptotic corrections. Clearly the asymptotic behaviour sets in only for very large 
values of p. Looking at the data with pmin = 3.0, one sees that the estimates show a 
systematic trend as A^ increases. If we analyze together all data with A^ > A^min we 
have (pmin = 3.0) 

-0.184 ±0.004 

-0.173 ±0.005 ..„.n-^ , .g^gs 



0",, 



-' 'mm 


= 1000 


A^ ■ 

^ 'mm 


= 2000 


N ■ 

-' 'mm 


= 4000 


A^ ■ 

^ 'mm 


= 8000 


N ■ — 

-' ' mm 


1000, 


N ■ — 

-' ' mm 


2000, 


N ■ — 

-' ' mm 


4000, 


N ■ — 

-' ' mm 


8000. 



-0.169 ±0.005 
-0.157 ±0.006 
0.289 ±0.002 
^ .0.284 ±0.002 .,^,n--www, , . 

Jw,oc 1 0.282 ±0.002 ^^ - /innn l^' ') 

0.278 ± 0.002 

The presence of confluent corrections and of non-asymptotic terms of opposite sign 
makes difficult to evaluate reliably a^ and fw,oo'- clearly large values of N are needed 
to see the scaling regime and large values of p are required to observe the correct 
asymptotic behaviour. From the results of the fits above we can conclude that there 
is a reasonable agreement between the theoretical estimates and the numerical results 
although a precise quantitative check is difficult. 

Let us now consider the distribution function /(p). In order to determine this 
function, we have computed the probability Pn{''i^) from the Monte Carlo data. A 
graph of this quantity as a function of r^ shows strong oscillations due to the un- 
derlying lattice structure. In order to reduce these effects, we have used a procedure 
analogous to that used in Ref. [0. Given a number Ngh, we define r^ = nNgh and 
an averaged distribution function 

PNirn) = ^ E ^^W, (5.18) 

where Nn is the number of lattice points in the shell r'^_i < r^ < r^. For Ngh fixed, 
in the scaling limit |r| — >■ cxd, A^ —>■ oo, with p fixed, P^{r) converges to /(p), so that 
one can use the distribution (|5.18| ) in order to compute the EEDF. The advantage is 
that lattice oscillations disappear in the averaging procedure. Of course, one should 
always check that the results do not depend on Ngh. As expected, as long as the 
number of points falling in each shell is sufficiently large and y/Ngh -C ^n-, the final 
estimates are not sensitive to A^s^. 

We have closely repeated the analysis performed for the wall-to-wall EEDF. The 
final results are in reasonable agreement with the theoretical predictions. First, we 
have performed two different sets of fits in order to determine D and 5. As before, 
we consider 

log(/(p)) = log/oo-/^p^ (5.19) 

log(p-°-2^^ /(p)) = log/oo-V, (5.20) 



for various values of pmin- The second fit keeps into account the presence of p" using 
the theoretical prediction for cr, Eq. ([4.39|) . The results of the fits for fixed values of 
A^ and different Ng^ are reported in Table ||. No significant dependence on iVmin and 
Pmin is visible in these results. Considering all data with A^ > 1000, and using for 
each A^ the largest Nsh appearing in Table ^, we obtain from the first fit 



D = 
while the second one gives 
6 = 
D = 



2.504 ±0.002 


Pmin = 3.0, 


(5.21) 


2.481 ±0.004 


Pmin ^ O.O, 


0.1222 ±0.0004 


Pmin = 3.0, 


(5.22) 


0.1277 ±0.0008 


Pmin ^ O.O, 


2.444 ± 0.002 


Pmin ^ 3.0, 


(5.23) 


2.441 ± 0.004 


Pmin ^ O.O, 


0.1397 ±0.0005 


Pmin ^ 3.0, 


(5.24) 


0.1406 ±0.0009 


Pmin ^ 3.5. 



Fit (|5.19|) gives estimates that show strong scaling corrections, clearly due to the 



neglected power term. The asymptotic values are difficult to estimate from this fit. 
In any case, we should observe that the estimates have the correct trend towards the 
expected results. The second fit is more stable. The estimates are in much better 
agreement with the theoretical results, although larger values of pmin are necessary 
to confirm the theory at the level of the statistical precision we have here. From the 
second fit we can also estimate f^o- We obtain 

^r 0.0158 ±0.0001 Pmin = 3.0, 

•^°° 10.0159 ±0.0001 pi„in = 3.5, ^^^ 



in agreement with Eq. ( [4.411) . 



Finally, in order to obtain estimates of /oo and a, we have performed a fit of the 
form 

log [/(p) exp(Dp^)] = log U + a logp, (5.26) 

using the theoretical estimates for D and 6, for various values of p > pmin (see Table 
H). These fits become rapidly unstable with increasing pmin- Nonetheless, as we shall 
see, the final results are in reasonable agreement with the theoretical estimates even 
if one considers 1 < pmin < 2. Fitting all data with N > 1000, we obtain 

0.2454 ± 0.0003 
a = ( 0.2426 ± 0.0004 
0.2351 ±0.0006 
0.01608 ±0.00001 
/oo = { 0.01612 ±0.00001 
0.01625 ±0.00002 

These results are not far from the theoretical estimates, although they show a trend 
with increasing pmin which is the opposite of what one expects theoretically: indeed 
the difference between the numerical and the theoretical estimates increases with 

19 



Pmin = 1-0, 




Pmin ^ 1-5, 


(5.27) 


Pmin = 2.0, 




Pmin = 1.0, 




Pmin = 1-5, 


(5.28) 


Pmin = 2.0. 





Pmin — >^ cxD. It shouM be noticed however that this behaviour could be a result of 
corrections to scaling: indeed, for pmin = 2, one observes that the estimates of o 
increase with A^. If we analyze together all data with A^ > A^min, we have (pmin = 2) 



I U.ZOUr IC U.UUUU l-^mm = ^UUU, , s 

^ ~ \ n OQ7n -I- n nnnv at . _ /innn (^O.zyj 



0.2351 ±0.0006 


iVmin = 


1000, 


0.2367 ±0.0006 


A^min = 


2000, 


0.2379 ±0.0007 


A^min = 


4000, 


0.2393 ±0.0009 


TV ■ — 

^ 'mm 


8000, 


0.01625 ±0.00001 


^ "mm 


= 1000, 


0.01622 ±0.00001 


-' »min 


= 2000, 


0.01619 ±0.00001 


A^min 


= 4000, 


0.01616 ±0.00002 


TV ■ 

^ 'mm 


= 8000. 



f = J "-"-^^^^ -^"-""""-^ ^'mm — ^""v., {^'\^\ 

•'°° 1 n nifiin -I- n i^[\i^i^^ at . — /innn y^-^^jj 



Therefore, we observe two opposite effects: a increases with increasing pmjn because 
of non-asymptotic corrections in /(p), while it decreases with increasing TV^jn because 
of corrections to scaling. Such a behaviour is not unexpected, since /^(p) was found 
to behave in exactly the same manner. For these reasons an accurate numerical check 
of the predictions for a and /oo is difficult: it is however reassuring that all estimates 
are reasonably near the theoretical results. 

Finally, we have studied the behaviour of /(p) for p ^ 0. In this case we have 
performed fits of the form 

log/(p) = log/i±^logp, (5.31) 

using, in each case, only data with pmin < p < Pmax- We have introduced here two 
cuts, Pmin and Pmax- The meaning of the latter is clear: it plays the role of pmin 
in the analysis of the large-p behaviour of the EEDF. The second parameter pmin is 
introduced to eliminate spurious lattice effects. Indeed, the scaling limit is obtained 
taking |r| to infinity. In other words, small values of r should be discarded. In our 
fits we have taken pmin = 0.1. Of course, this choice introduces a bias, and one should 
study the hmit pmin — > to obtain the correct asymptotic behaviour. In our case, the 
systematic error appears to be small: indeed the estimates are stable with respect 
to small changes of this parameter. The results for fixed values of TV are reported in 
Table 0. The data show a small systematic variation with TV. Using all data with 
TV > 2000, we obtain 



0.223 ±0.003 


Pmax = 0.6, 






0.253 ±0.006 


Pmax U.4, 




(5.32) 


0.281 ±0.031 


Pmax ^ U.Z, 






0.01539 ± 0.00004 


Pmax 


0.6, 




0.01599 ±0.00012 


Pmax 


0.4, 


(5.33) 


0.01680 ±0.00096 


Pmax 


0.2. 





/l 



The estimates of 9 show a systematic variation with Pmax, indicating the presence 
of strong non-asymptotic corrections. For pmax ~ 0.40, 9 is in reasonable agreement 
with the theoretical prediction. The constant /i shows a similar trend, approaching 
the theoretical value for pmax ~ 0.40. From /i we can compute /o = /i2^. Using 
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the theoretical prediction for 6, we have /o = 0.01853(5), 0.01925(14), 0.0202(11), 
corresponding respectively to pmax = 0.6, 0.4, 0.2. These estimates are in reasonable 
agreement with the theoretical results presented before that predicted 0.015 ^ /o ^ 
0.019. 

In conclusion, our Monte Carlo results confirm the theoretical results of the previ- 
ous Sections. Notice that the theoretical predictions are more precise than the Monte 
Carlo estimates, in spite of the large statistics and of the very long walks used. 
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A Definitions and properties of the basic functions 

In this Appendix we report the definitions and asymptotic expansions of the functions 
that appear in our two-loop computation of the end-to-end distribution function. 

Function Ni{x) 

We define _ 

A^i(a;) = / dt(l-xt) e~^'* logt log|l-t|. (A.l) 

Jo 

The expansion of Ni{x) for large values of x is easily computed. Indeed, in this limit 
the relevant contribution is due to the region t ^ 0. It is then enough to expand 
log |1 — t| in powers of t and integrate term by term. One obtains 

^i(^) = E ^—^ H(^ + 1) + 1] - logo: E -^^, (A.2) 

n=l -^ n=l -^ 

where iplx) is the logarithmic derivative of Euler's F-function. To obtain the asymp- 
totic expansion for small values of x, one first notices that A'^i(x) satisfies the differ- 
ential equation 

N[{x) + Ni{x) = ^\2--fE-logx- e-"Ei(x)l , (A.3) 

X^ L J 

where 7^; ~ 0.5772156649 is Euler's constant, and Ei is the exponential integral func- 
tion ||50|. Solving the previous equation, one obtains a different integral representation 
for Ni{x): 



r^ (if r 

mix) = e-M - e* {2-jE- logt) - Ei(t) - (2 + t)(l - 7^ - logt) 

2 
(7ij + logx) + (1 - 7^) logx - i log^ x + l + 'jE- ^TT^ - ^7! 
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• (A.4) 



Using this expression it is trivial to obtain the expansion of A'^i(x) for x —* 0: 



iVi(x) 



-(logo; + 7i5) + i [2 - 27s - 7I - Ivr^ - 2(1 + -/e) logx - log'x 



X 
X 

+ 2 



I - IE + ll + 7-7^^ + {2lE - l)'^ogx + log^ x] +0(x2logx). (A.5) 



Function N2{x) 
We define 



N2{x) 



dt 



dz 2;3/2ga;2 



t'^^^^hc2m{z + lY 



Kl[^t]-HJ^t 



(A.6) 



where C is a loop contour going counterclockwise around the negative real axis, Ki (z) 



and Ji{x) are Bessel functions [0, and H^if) is defined by 



Hr{t) 



-rt 



+ r^ + 37i5 + 3 log ■ 



t3 + t2 + 2t (v 2 ' ^ 
3r r^ 3r 3r t 



t 



(A.7) 



Note that Hrif) is such that Klit) - Hr{t) ~ tlog^t for |t| -^ 0. 

We want now to derive the asymptotic behaviour of N2{x) for large and small 
values of x. Substituting w = xz and s = tj ^fx we can rewrite 



°° ds 



JO s Jc 2m (w + x^ 



du; tiJ^/^e"' 



ir-' (v^s) - //, (Vws)] . (A.8 



To derive the small-x behaviour, we expand Ji{sy/x) (the corresponding series con- 
verges everywhere since Ji{z) is an entire function) obtaining 



N2{x) = E 



(-x)" f^dsfsY'^+'^ f dw w^/2 



^on!(n + l)!io sV2/ Jc2-Ki{w + xY 



Kl (v/^s) - iT, (v^s) 



(A.9) 

Now, Ki{z) ~ 2; ^/^e ^ for |z| ^ 00. Since |argw| < tt, we can define t = ^/ws, and 
rotate the contour so that t belongs to the positive real axis. We obtain 



N2ix) = y: 



i-xy 



^ dt ft 



2n+l 



,^on!(n + 1)170 t \2 
The last integral can be done exactly. We get 
dw ^^-"e"' 



K\{t) - Hrit) 



dw w^ "e*" 
c 2m {w + xY 



(A.IO) 



c 2ni (w + xY 



n-2 

X '"e "[x + n — 1) + y^ 

k=0 



1 \n+l^—n^—x 



n — k — 1 



-X 



\k~n 



(A.ll) 
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with the convention that, for n < 2, the summation is zero. Substituting in Eq. 
( lA.lOp we obtain 



N2{x) 



roo fjf V f 

-e-^'l j[KKt) - Hr{t)] [-/o(t) + (x - 2)h{t) 



+ E- 

A;=0 



-X) 



k\ 



°° dt 
T 



[Kl{t) - H^{t)] 



X 



+ ^ n — h — A / f \ 2n+l 



(A.12) 



This expression can be simphfied and one obtains the final result 



N,{x) = y: 



k=0 



k\ 



E 



"on!(n + l)! 



, POO fif / + \ ^n-ti 

-j^ j[Kl{t)-HAt)]{^^ . (A.13) 



This expansion converges absolutely for all values of x and allows the computation 
of N2{x) up to quite large values of x with a small effort. It gives immediately the 
small-x expansion of N2{x). 

To compute the large- a; behaviour of N2{x), we start by introducing the Mellin 
transforms of [Kf{t) — Hr{t)]/t and of Ji{t). Explicitly, we define 



M{z)= / dtt'-^ Kt(t)~Hr{t) 

Jo L 



and compute 



^ ^+1 ' 






(A.14) 



(A.15) 



Eqs. ( |A.14|) and ( [A.15|) are defined for Re^; > and — 1 < Rez < | respectively. 
Using the Parseval formula for Mellin transforms, we can rewrite N2{x) as 



N2{x) = x'^l^ 



c 2'ni[w -\- xV J-i~ioo 2711 r 0^ 



where 7 is any real number with — 1 < 7 < 1. Let us now show that N2{x) 
a; — >■ 00 faster than x~^~^'^, for any e > 0. Indeed 



where 



|iV2(x)| < a;-^/2-^/22^-^ 



dy 



\dw\ 



,Reui 



c 2tt \w + xP 



\w 



(7+3)/2 



/H, 



/H= r-^|M(l-7-^l/)|e- 

J-oo ZTT 



y{a.Tgw)/2 



-m+w+^\ 


^y 2 1 


p/3- 


~i-w\ 


j- 1 


2 > 



(A.16) 
-^ for 

(A.17) 
(A.18) 



Let us first consider the integral f{w). For \y\ — > 00, an easy computation shows that 
|M(1 — J — iy)\ ~ \y\P^'^^e~'^y^'^, where ^(7) is an exponent we do not need to know 
explicitly. In the same limit the ratio of F-functions behaves as |i/|'''~^. Thus, the 
integral f{w) is finite for |argw| < tt. For |argu7| -^ vr, depending on the value of 7, 
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the integral may be finite or diverge as a power of (|arg w| — vr). But, witli our clioice 
of C, |argw| — >■ TT corresponds to Rew — ^ — oo. In this hmit the integrand in ( [A.17| ) 
decreases exponentially, therefore making the whole expression finite. Finally notice 
that 



\w + x\ 



4 
< — 



X 



3x ^ X 

< Rew < — 

2 - - 2 



+ 



1 



(Imw) 



-,x 



3a; ^ X 

— < Rew < — 

2 - - 2 



where x (condition) is 1 if the condition is satisfied, otherwise. It follows 

|iV2(a;)| < x-5/2-7/227+1 ( \M 

+X-1/2-7/227- 



(A.19) 



c 27r 
1 / \dw\ 



.gRe«,i |(7+3)/2 



^Reio 



c 27r ilvnwY 



\w\ 






(A.20) 



where C is the part of C with — ^ < Rew < — |. It is possible to choose C so 
that |Imu7| is constant along C. It is then trivial to show that the second integral in 
Eq. (|A.2CI|) is bounded by x^e~^/^, where p is an appropriate power. Thus |iV2(a:)| < 
const x~^l'^~^l'^. Since 7 is arbitrary with 7 < 1, the result follows immediately. 
A little more work, using the technique presented for Nj^ix)^ allows to show that 
^"2(3^) ~ 0(2;~^log^a;). It is interesting to notice that, by adding additional terms in 
Hr{x), one can make N2{x) decrease faster: if Hr{x) is such that Kf{x) — Hj.{x) ~ 
0{x'^ log x), then N2{x) decreases faster than x~"/^"^"*"^, for any e > (more precisely 
N2ix) ~0(x-"/2-2iog3a;)). 

Function Ns{x) 
We define 



Nsix) 



dt 
T 



Ji{t)-¥+Tf 



'' ""''^E.iy-z. 



c 2m [z + 



zt 



(A.21) 



where C is a loop contour going counterclockwise around the negative real axis, Ji (x) 
is a Bessel function [Q and Hr(t) is defined in Eq. ( [A.?] ). 

The small-x behaviour is easily computed using the same procedure as before. We 
obtain 



Nsix) = y: 



-X 



^^ ^ k + l-n 



k\ 



E 



-2n-l 



k=2 - n=2^K^+l)!' 

where Hj.{z) is the Mellin transform of Hr{t)/t: 



Hr{2n + 2) 



(A.22) 



HJz) 





2-z 



dt^-'Hrit) 
r'^V{z - 4) + r'^V{z - 3) + \V{z - 2)(-3 + 67^ + 2r'^ 



+ ^T{z - l)(-9 + l^^E + 2r^) + \T{z - 2){i,{z - 2) - log2r) 



+ lT{z-l){i){z-l)-\og2r) 



(A.23) 
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This expansion converges absolutely for all values of x and it gives immediately the 
small-x expansion of N^lx). 

Let us now compute the behaviour for large values of x. We rewrite N^lx) as 



Ns{x) =Li(x)+L2(x), 



(A.24) 



where 



Li(x) 



X 



1/2 r ^ [j,{^s) + 6(1 - V^S) (-isxV2 + _Ls3^3/2 

r°° ds 



X 



dw vs^^'^e^ 



3^3/2 



X 



c 2m (w + x)2 



Hr{y/ws), 



(A.25) 



(A.26) 



and 9{x) is Heaviside's step function. 

Let us first consider Li{x). Using the Parseval formula for Mellin transforms we 
obtain 



Li{x) = X 



-1/2 



dw w^/^e"' /■7+ioo dz 



c 2711 {w + xy Jj~ioo 27ri 



w'/^x-'/^H{l - z)B{z), (A.27) 



where 



Biz) = r-'J^ - h^ + ' 



r(^) 2^+1 16^ + 3' 



(A.28) 



and 7 is a real number satisfying —5 < 7 < —3. Now, rewrite the previous expression 
for Li as 



Li(x) 



X 



-1/2 



dw w^/^e"' n+ioo d^ 



c 2m [w + xY J-y-i 



J'y—ioo Zm 



—X 



_i/2 / an^ w e ^ Res\w^/^x'^/^H{l~z)B{z)\, (A.29) 



c 2m {w + xy 



n=—3 



where, in the first integral, 7 is real such that < 7 < 1. Repeating the discussion 
presented for N2{x) one can easily show that the first integral in Eq. ([A.29| ) behaves 
as 0{x~^~^'^). Neglecting terms of this order, we obtain 



L,W = _(-41og2 + 4,, + ll)/^-^-^-j, 

27 



128 



27 3 2 -^7, 3,2 ^ 1 \ /■ 

+ 50. -5 log 2) y^ 



dw we^ 
c 2m {w + xy 



c 2m {w + x)'' 



log(w/x). 



(A.30) 
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Using Eq. ([7011) and (n > 0) 



dw w^'e'" , (-l)"+^n! / 2(n + 11 
' ;^^^-^logw;= ^ ^ (1 + ^ ^ 

C iTTi [W + XY 



0(x- 



x^ 



X 



we obtain 



Liix) 



27 3 



-o7s+olog2 — + 0(x Mogx 



(A.31) 



(A.32) 



,128 8 8 " / x^ 

The computation of the asymptotic behaviour of L2{x) is completely analogous. We 
have, discarding terms of order x~^ log^ x, 



r . N -1/2 /■ ^^ ^'^'e'^ ^ T> 

L2\x) = —X ' / rr > Kes 

^ ^ ic 27ri (w + x)\t^3 ^=" 

Using the previous results and (n > 0) 



^./2^-./2^^^ - Z) 



11 11 



2z + l 16^ + 3 

(A.33) 



dw w^e'" 2 _ 2(-l)"+V("'+l)^! ^/ -3x 

C 2TTi (w + xY x'^ 



we obtain 

L2{x) 



32x2 
Summing up, we have 

Nsix) 



X H h 67E — 12 log 2 — 6 log X 



32x^ 



(A.34) 



+ 0(x-'^log^x). (A.35) 



[x + 8-67s-61ogx] +0(x"''log^x). (A.36) 



Fourier transform of N2{Q'^) + N^IQ^) 

We want now to report the small-p behaviour of the Fourier transform of A^2(<5^) + 
N^IQ"^) that we use in Sec. ^. We define 



N{p)^ J^^e^^-^{N2{Q') + N,{Q^^ 



(2vr) 



where, cf. Eq. ( |3.15| ), 

iV2(x) = Ar2(x) - e-"(l + x)Ar2(0) - e-"Ar'2(0). 
Using the definitions of N2{x) and A^3(x) we can rewrite 



(A.37) 



(A.38) 



Ic 27ii {w + x)2 Jo s 
Using the fact that, for any function h{Q'^), we have 



Ar2(x)+Ar3(x) = x'^'^ ^ ^^^^^^, C -KK^V^) IMV^s) - ^V^s + ^^x'^h' 

(A.39) 

(A.40) 



r d^O 1 /"oo 



Att^P Jo 
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we can perform the integral over Q, cf. formula 6.541 of Ref. |E0[, obtaining 



N{p) 



IGtt^P Jc 27ri 



dw 



-.we 



ds 



s 



Kl{s) 



K^{V^p)(l2{s) + h{s)-l-ls' 



+ [K^{V^p) + K,{^p)) (his) -Is- ^s^ 



+ttV / P^^^"" 1-pmpv^) + hipv^)] r -Kt 

loTT'^p Jc 2m [ Jp s 

-p[K,{p^) + K,{p^)] / -Kl{s^)h{s^\ 

J p s 

/•oo 

+Ki{py/^) / dsKf{sVw)[l2{sVw) + loisVw)] 
Jp 

+/i(pv^) / dsKUsV^)[K2is^/^) + Kois^/^)] 
Jp 

Expanding for p — * 0, after a lengthy calculation, we obtain 

1 3.0-0 1 



'(.v^) 



(A.41) 



N{p) 



■^°S 2" 128vr^ 



1287rV 1287r2^"^ ^ i 98^2 v^ + ^Ts) log 



7b 



1024 2567r2 2567r2 5127r2 



7i; + 



A 



327r2 



0(p2 logV),(A.42) 



where 



A= r ds Kl{s){l + \s^) -e 



s/2 



1 1 3 3 3 , s 

+ ^r^ + —lE - — + — lo. 



2s2 2s 8s 2s 2 

Numerically A ^ -1.30204. 
Function gi{x) 
We define 

gi{x) = i(l - e"^) + 1(1 - x)e-^ (Ei(x) - logx) , 



(A.43) 



(A.44) 



where Ei(x) is the exponential integral function [^. We are interested in the be- 
haviour of gi{x) for large and small values of x. These expansions are easily obtained 
using the corresponding results for Ei(x). For x ^ 1 we obtain 



^ n=0 "-• 



For X ^ 1 we have the asymptotic expansion 



1 ^ (n + l)(n + l)! 

^i(^) =-92^ :;:^ • 



(A.45) 



(A.46) 
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Function g2{x) 
We define 



^^^"^ = Tx 



e ^Ei(a;) + 7^; + logx log a; 



+ -e ""(1 -x) logx [logx - 2Ei(x)] 



1 



21oga;-4Ei(a;) -7r2(l-a;) - -(1 - a;)A^i(a;). (A.47) 



For a; <^ 1 we have 
92{x) = 

For a; ^ 1 we have 

Function g^{x) 
We define 

For a; ^ 1 we have 



oE(^ + i) 



V'2(n+ 1) - V^'(n + 1) 



-X] 



n=0 



n! 



\ (^^ - i-') + I (-2 + 27i. - 7l + Ivr^) + 0(a:^). 



, , 1 ^ A;!A;^(A; + 1) 
^2W = tZ^ 



k=l 



X 



k+l 



9z{x) = —1+xe ^Ei(a;) + e ""(1 — a;loga;). 



93[X) 






X 



For X ^ 1 we obtain 
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Figure 1: Graphs contributing to the EEDF to order e^. 
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Figure 2: Plot of /(p) for different values of A^. 
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1.51169 ±0.00043 
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2.9833 ±0.0019 
2.9933 ±0.0020 
3.0011 ±0.0015 
3.0053 ±0.0026 
3.0158 ±0.0042 


7.044 ±0.011 
7.1600 ±0.0098 
7.2360 ±0.0096 
7.2811 ±0.0095 
7.3153 ±0.0075 
7.334 ±0.012 
7.386 ±0.023 
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20.306 ±0.050 
20.647 ±0.050 
20.849 ±0.050 
20.996 ±0.039 
21.091 ±0.066 
21.34±0.11 


M* 
B 
A 

CL 


1.51406 ±0.00089 

-(0.00314 ±0.00084) 

0.538 ±0.079 

70% 


3.018 ±0.005 

-(0.016 ±0.0046) 

0.555 ±0.085 

66% 


7.387 ±0.023 

-(0.070 ±0.021) 

0.572 ±0.092 

67% 


21.34±0.13 
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th 
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7.392 ±0.015 


21.35 ±0.06 


ph 


1.50876 ±0.00023 


2.993 ±0.001 


7.292 ±0.005 


20.94 ±0.02 



Table 1: Monte Carlo results for the invariant ratios M2k- "CL" is the confidence 
level of the fit Msfc = M^ + S2fc(A^/8000)-^. "th" is the theoretical prediction ( ^3TD 
and "ph" the phenomenological prediction ( [4.52|) . 
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Table 2: Monte Carlo results for the invariant ratio R2k- 
prediction. 
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2.479(13) 
0.1318(30) 
0.2747(50) 


2.423(30) 
0.1462(81) 
0.305(17) 


2.444(74) 
0.141(20) 
0.293(50) 


2000 


s 

D 

Jw,oo 


2.420(13) 
0.1483(34) 
0.2506(48) 


2.434(29) 
0.1442(78) 
0.244(13) 


2.476(84) 
0.133(21) 
0.224(43) 


2.463(13) 
0.1346(31) 
0.2735(50) 


2.468(29) 
0.1336(72) 
0.272(14) 


2.494(83) 
0.127(20) 
0.259(48) 


4000 


s 

D 

J 10,00 


2.407(13) 
0.1517(33) 
0.2541(47) 


2.447(29) 
0.1408(75) 
0.235(12) 


2.456(74) 
0.138(20) 
0.232(40) 


2.452(13) 
0.1375(31) 
0.2767(50) 


2.483(28) 
0.1297(68) 
0.261(13) 


2.486(75) 
0.129(19) 
0.262(44) 


8000 


5 
D 

Jw,oo 


2.4139(98) 
0.1490(25) 
0.2464(35) 


2.419(23) 
0.1477(61) 
0.244(10) 


2.450(58) 
0.139(15) 
0.224(30) 


2.4590(99) 
0.1350(23) 
0.2684(37) 


2.454(23) 
0.1362(57) 
0.271(11) 


2.480(59) 
0.129(15) 
0.253(33) 


16000 


s 

D 

Jw,oo 


2.370(21) 
0.1617(61) 
0.2658(86) 


2.423(51) 
0.147(14) 
0.241(22) 


2.36(15) 
0.167(48) 
0.29(10) 


2.414(21) 
0.1465(55) 
0.2883(88) 


2.458(51) 
0.136(13) 
0.268(24) 


2.37(15) 
0.160(45) 
0.33(11) 


32000 


s 

D 

Jw,QO 


2.433(32) 
0.1435(79) 
0.236(10) 


2.382(83) 
0.157(24) 
0.257(39) 


1.69(22) 
0.62(28) 
2.4 (2.3) 


2.482(33) 
0.1292(73) 
0.257(11) 


2.422(83) 
0.144(22) 
0.282(41) 


1.74(22) 
0.56(25) 
2.3 (2.1) 



Table 3: Results for the fit g{p) = /^^ooexp(— Dp^) for p > pmin- 
columns refer to g{p) = fw,Mc{p), the last three columns to g{p) - 
Here U,mc{p) is the Monte Carlo "wall-to-wall" EEDF. 



The first three 



N 




Pmin = 2. 


Pmin — 2.5 


Pmin — 'J- 


500 


CTw 
Jw.oc 


-0.1153(31) 
0.27406(81) 


-0.1752(63) 
0.2934(20) 


-0.252(12) 
0.3232(50) 


1000 


O'w 
Jw.oc 


-0.12149(24) 
0.27204(62) 


-0.1641(47) 
0.2855(15) 


-0.234(10) 
0.3117(38) 


2000 


O'w 
Jw.oc 


-0.1254(24) 
0.27076(60) 


-0.1527(45) 
0.2792(14) 


-0.1902(97) 
0.2926(35) 


4000 


Ow 
Jw.oc 


-0.1268(22) 
0.26937(56) 


-0.1513(44) 
0.2770(13) 


-0.1954(94) 
0.2927(34) 


8000 


Ow 
Jw,oo 


-0.1280(17) 
0.26850(44) 


-0.1417(34) 
0.2727(10) 


-0.1537(72) 
0.2768(24) 


16000 


(Jw 
Jw.oc 


-0.1314(34) 
0.26874(87) 


-0.1473(70) 
0.2735(21) 


-0.2001(17) 
0.2919(60) 


32000 


Ow 
Jw.oc 


-0.1087(51) 
0.2624(12) 


-0.125(11) 
0.2670(33) 


-0.090(33) 
0.256(10) 



Table 4: Results for the fit fn,,Mc{p) exp{Dp^) = /^,ooP'^"' for p > Pmin- fw,Mc{p) is 
the Monte Carlo "wall-to-wall" EEDF. D and 5 have been set equal to the theoretical 
predictions, Eqs. (^) and {^^^. 
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N 


Nsh 




/mc(p) 


fMcip) ■ P-'--'"' 








Pinin ^ O 


Pmin — O.O 


Pmin — 4 


Pmin =^ 'J 


Pmin — 3.5 


Pmin = 4 


500 


1 
10 


s 

D 

/oc 

s 

D 

/oc 


2.5368(59) 
0.1170(12) 
0.01827(15) 
2.5662(59) 
0.1113(12) 
0.01762(14) 


2.479(11) 
0.1305(27) 
0.02045(43) 

2.545(11) 
0.1158(24) 
0.01833(36) 


2.296(23) 
0.1877(86) 
0.0333(21) 
2.464(23) 
0.1359(61) 
0.0224(12) 


2.4733(58) 
0.1346(14) 
0.01583(14) 
2.5019(58) 
0.1282(13) 
0.01523(13) 


2.430(11) 
0.1461(30) 
0.01732(38) 

2.495(11) 
0.1299(26) 
0.01545(32) 


2.260(23) 
0.2050(93) 
0.0279(18) 
2.425(23) 
0.1490(66) 
0.0185(11) 


1000 


1 
20 


6 
D 

/oc 

5 
D 

/oc 


2.4892(47) 
0.1262(10) 
0.01905(13) 
2.5336(47) 
0.11696(96) 
0.01800(12) 


2.4078(88) 
0.1473(24) 
0.02244(39) 
2.5094(88) 
0.1224(20) 
0.01885(30) 


2.202(18) 
0.2216(81) 
0.0395(21) 

2.461(19) 
0.1345(51) 
0.02107(97) 


2.4268(46) 
0.1450(12) 
0.01655(12) 
2.4698(46) 
0.1346(11) 
0.01560(11) 


2.3601(89) 
0.1646(27) 
0.01911(35) 
2.4596(87) 
0.1373(21) 
0.01593(27) 


2.167(18) 
0.2419(88) 
0.0335(19) 

2.422(18) 
0.1477(52) 
0.01744(79) 


2000 


40 
80 


s 

D 

/oc 

D 

/oc 


2.5123(44) 
0.12072(94) 
0.01822(12) 
2.5161(44) 
0.11993(94) 
0.01813(11) 


2.4940(83) 
0.1249(19) 
0.01885(29) 
2.5035(83) 
0.1228(19) 
0.01856(28) 


2.421(17) 
0.1443(48) 
0.02257(95) 

2.445(17) 
0.1376(46) 
0.02135(89) 


2.4489(43) 
0.1389(11) 
0.01580(10) 
2.4526(43) 
0.1380(11) 
0.01573(10) 


2.4442(81) 
0.1401(21) 
0.01595(25) 
2.4534(83) 
0.1378(21) 
0.01569(25) 


2.382(17) 
0.1585(53) 
0.01875(82) 

2.406(17) 
0.1513(50) 
0.01771(76) 


4000 


100 
200 


5 
D 

/oc 

D 

/oc 


2.5031(42) 
0.12229(91) 
0.01828(11) 
2.5071(43) 
0.12145(92) 
0.01818(11) 


2.4716(79) 
0.1297(19) 
0.01940(29) 
2.4823(79) 
0.1272(19) 
0.01905(28) 


2.4289(16) 
0.1413(45) 
0.02168(85) 
2.457(16) 
0.1338(43) 
0.02036(79) 


2.4398(41) 
0.1407(10) 
0.01587(10) 
2.4439(42) 
0.1397(10) 
0.01578(10) 


2.4221(79) 
0.1455(22) 
0.01644(26) 
2.4325(78) 
0.1427(21) 
0.01613(25) 


2.390(17) 
0.1552(52) 
0.01799(78) 

2.418(17) 
0.1471(49) 
0.01685(72) 


8000 


200 
800 


(5 
D 

/oc 

D 

/oc 


2.4953(35) 
0.12364(78) 
0.018323(94) 

2.5003(36) 
0.12258(78) 
0.018205(94) 


2.4839(67) 
0.1263(16) 
0.01872(23) 
2.4964(68) 
0.1235(16) 
0.01833(23) 


2.447(14) 
0.1358(36) 
0.02048(67) 

2.484(14) 
0.1266(34) 
0.01890(60) 


2.4320(35) 
0.14230(88) 
0.015919(86) 

2.4369(35) 
0.14108(88) 
0.015808(87) 


2.4341(65) 
0.1418(17) 
0.01585(20) 
2.4462(66) 
0.1386(17) 
0.01551(20) 


2.408(14) 
0.1493(42) 
0.01696(61) 

2.444(14) 
0.1394(37) 
0.01561(52) 


16000 


800 
1200 


(5 
D 

/oc 

5 
D 

/oc 


2.4913(53) 
0.1245(12) 
0.01842(14) 
2.4935(53) 
0.1240(12) 
0.01837(14) 


2.4619(99) 
0.1316(24) 
0.01950(36) 
2.466(10) 
0.1305(25) 
0.01936(37) 


2.409(22) 
0.1461(62) 
0.0222(12) 
2.423(22) 
0.1421(60) 
0.0215(11) 


2.4283(52) 
0.1432(13) 
0.01601(13) 
2.4306(52) 
0.1426(13) 
0.01595(13) 


2.4125(97) 
0.1475(27) 
0.01654(32) 
2.4169(98) 
0.1464(27) 
0.01641(32) 


2.370(20) 
0.1605(63) 
0.01847(96) 

2.384(20) 
0.1562(62) 
0.01786(92) 


32000 


5000 


(5 
D 

/oc 


2.4837(87) 
0.1258(19) 
0.01844(23) 


2.479(16) 
0.1268(38) 
0.01855(55) 


2.366(33) 
0.159(10) 
0.0246(21) 


2.4208(85) 
0.1446(22) 
0.01603(21) 


2.429(16) 
0.1422(42) 
0.01570(49) 


2.327(32) 
0.174(11) 
0.0205(18) 



Table 5: Results for the fit g{p) = foo exp{—Dp^) for p > Pmm- 
refer to g{p) = fMc{p), the last three columns to g{p) = /mc 
is the Monte Carlo EEDF. 



The first three columns 
(p)p-o-255. Here /mc(p) 
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N 


Nsh 




Pinin -L- 


Pmin — 1.0 


Pniin — ^- 


500 


1 
10 


a 

loo 

a 

foo 


0.25638(80) 
0.015976(11) 
0.25607(81) 
0.015976(11) 


0.2468(11) 
0.016130(17) 

0.2463(11) 
0.016132(17) 


0.2237(17) 
0.016540(29) 

0.2230(17) 
0.016548(29) 


1000 


1 
20 


a 

foo 

a 

foo 


0.25146(64) 
0.016019(9) 
0.25030(65) 
0.016027(9) 


0.24592(89) 
0.016108(13) 
0.24418(90) 
0.016126(14) 


0.2304(14) 
0.016382(23) 

0.2276(14) 
0.016420(23) 


2000 


40 
80 


a 

foo 

a 

foo 


0.24691(61) 
0.016062(9) 
0.24688(63) 
0.016062(9) 


0.24285(85) 
0.016128(13) 
0.24276(87) 
0.016129(13) 


0.2324(13) 
0.016311(22) 

0.2325(13) 
0.016310(22) 


4000 


100 
200 


a 

foo 

a 

foo 


0.24530(59) 
0.016074(8) 
0.24532(60) 
0.016073(8) 


0.24121(82) 
0.016141(12) 
0.24115(84) 
0.016141(13) 


0.2351(13) 
0.016247(21) 

0.2350(13) 
0.016248(22) 


8000 


200 
800 


a 

foo 

a 

foo 


0.24274(50) 
0.016103(7) 
0.24278(53) 
0.016102(7) 


0.24240(69) 
0.016109(10) 
0.24238(73) 
0.016109(11) 


0.2376(11) 
0.016194(18) 

0.2375(11) 
0.016194(19) 


16000 


800 
1200 


a 

foo 

a 

foo 


0.24322(75) 
0.016092(11) 
0.24330(77) 
0.016090(11) 


0.2434(11) 
0.016088(16) 

0.2434(11) 
0.016088(16) 


0.2419(16) 
0.016116(27) 

0.2415(16) 
0.016122(27) 


32000 


5000 


a 

foo 


0.2413(13) 
0.016115(18) 


0.2412(18) 
0.016115(27) 


0.2440(28) 
0.016066(47) 



Table 6: Results for the fit /mc(p) exp(Dp^) = fooP" for p > pmin- fMcip) is the 
Monte Carlo EEDF. D and 5 have been set equal to the theoretical predictions, Eqs. 



(2.9) and ( 4.32 ). Data with p > 4.5 have not been included in the fit. 
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N 


Nsh 




Pmax = 0.8 


Pmax = 0.6 


Pmax = 0.4 


Pmax = 0.2 


500 


1 


9 

h 


0.2177(53) 
0.01470(5) 


0.2680(87) 
0.01542(11) 


0.303(19) 
0.01610(38) 


0.50(13) 
0.0228(54) 


1000 


1 


9 

h 


0.1975(41) 
0.01474(4) 


0.2453(67) 
0.01543(9) 


0.274(14) 
0.01600(27) 


0.339(55) 
0.0179(18) 


2000 


1 


9 

h 


0.1906(39) 
0.01479(3) 


0.2315(63) 
0.01539(8) 


0.270(13) 
0.01615(27) 


0.194(73) 
0.0140(19) 


4000 


1 


9 
h 


0.1895(37) 
0.01479(3) 


0.2258(60) 
0.01536(8) 


0.284(13) 
0.01654(26) 


0.343(61) 
0.0185(21) 


8000 


1 


9 

h 


0.1874(31) 
0.01489(3) 


0.2275(51) 
0.01546(7) 


0.251(11) 
0.01593(21) 


0.313(53) 
0.0177(17) 


16000 


1 

10 

20 


9 

h 

9 

h 

9 

h 


0.1721(46) 
0.01487(4) 
0.1745(46) 
0.01486(4) 
0.1750(46) 
0.01486(4) 


0.2000(74) 
0.01527(10) 
0.2035(75) 
0.01528(10) 
0.2045(75) 
0.01529(10) 


0.186(15) 
0.01497(29) 

0.191(16) 
0.01500(29) 

0.194(15) 
0.01506(29) 


0.182(75) 
0.0151(21) 
0.189(76) 
0.0152(21) 
0.199(77) 
0.0154(22) 


32000 


1 

10 

20 


9 

h 

9 

h 
9 

h 


0.1416(75) 
0.01482(7) 
0.1579(75) 
0.01480(7) 
0.1597(75) 
0.01480(7) 


0.178(12) 
0.01535(16) 

0.200(12) 
0.01542(16) 

0.203(12) 
0.01543(16) 


0.212(26) 
0.01609(51) 

0.243(26) 
0.01637(52) 

0.247(26) 
0.01640(52) 


0.28(13) 
0.0179(42) 

0.36(13) 
0.0198(47) 

0.38(13) 
0.0203(48) 



Table 7: Results for the fit /mc(p) 
EEDF. 



fip'^ for p < Pmax- fncip) is the Monte Carlo 
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